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Abstract 

We derive a new acoustic black hole metric from the Abelian Higgs model. In the non-relativistic 
limit, while the Abelian Higgs model becomes the Ginzburg-Landau model, the metric reduces to 
an ordinary Unruh type. We investigate the possibility of using (type I and II) superconductors as 
the acoustic black holes. We propose to realize experimental acoustic black holes by using spiral 
vortices solutions from the Navier-stokes equation in the non-relativistic classical fluids. 
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1 Introduction 

Black hole is characterized by the event horizon, a spherical boundary out of which even the light 
cannot escape. In 1974, Hawking announced that black holes are not black at all, and emit thermal 
radiation at a temperature proportional to the horizon surface gravity. Although several decades passed, 
the experimental verification of Hawking radiation is still elusive. 

Compared with the difficulties on the astrophysical side, analog models of general relativity is 
more appealing since these models are shedding light on possible experimental verifications of the 
evaporation of black holes. In the remarkable paper of Unruh [HI, the idea of using hydrodynamical 
flows as analog systems to mimic a few properties of black hole physics was proposed. In this model, 
sound waves rather than light waves, cannot escape from the horizon and therefore it is named "acoustic 
(sonic) black hole". A moving fluid with speed exceeding the local sound velocity through a spherical 
surface could in principle form an acoustic black hole. The event horizon is located on the boundary 
between subsonic and supersonic flow regions. 

In general, any fluid flows could be the candidate of acoustic black holes. Many classical and 
quantum systems have been investigated as black hole analogues, including gravity wave flSl, water 
BH, slow light [|5]-|7l, an optical fiber (jSJ, and an electromagnetic waveguide flU. But considered the 
detection of "Hawking radiation" of acoustic black holes, quantum fluid systems, rather than classical 
fluid systems are specifically preferred in that Hawking radiation is actually a quantum phenomenon. 
Up to now, superfluid helium II ifTOl . atomic Bose-Einstein condensates ifTTl . one-dimensional Fermi- 
degenerate noninteracting gas [[T2l have been proposed to created an acoustic black hole geometry in 
the laboratory. The possible application of acoustic black holes to quantum information theory has been 
suggested in [fT3l . 
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The first experimental realization of acoustic black hole reported was conducted in a Bose-Einstein 
condensate [[T4l . So far it has been assumed that the system is non-relativistic, because all the atomic 
system is likely in such regime. However, some strongly interacting system shows some evidences 
for exotic behavior that the matter is based on the heavy atoms but the basic excitations has massless 
character I115[|16L In fact the linear and massless dispersion relation is canonical for fermionic system. 
Therefore it is interesting to ask what happen to the acoustic black hole for the relativistic fluid system. 

In this paper, we explore a relativistic version of acoustic black holes from the Abelian Higgs model. 
Another motivation comes from the feasibility of detecting Hawking temperature associated with the 
acoustic black hole. Since the Hawking temperature depends on the gradients of flow speed at the 
horizon, detecting thermal phonons radiating from the horizons is very difficult. In fact the Hawking 
temperature calculated from models in Bose-Einstein condensates so far is very low (~ nano Kelvin). 
In constrast, the acoustic black holes created in the color superconducting phase in the dense quark 
matter and the pion superfluid phase at finite isospin density may provide us examples of relativistic 
superfluid and higher temperature acoustic black hole. The line elements of the black hole metric we 
obtained are found to be different from the non-relativistic version, which has been widely used for 
almost 30 years. 

2 Acoustic black holes from Abelian Higgs model 

The Abelian Higgs model is the Mexican-hat model coupled to electromagnetism with the action 

S = J d^xi-^F'^''F^,+ \{d^-teA^)<P\^ + m'\^\'-b\<P\^), (2.1) 

where F^,, = d^A^ — d^Ay. The Planck units fi = c = /cB = lare used here. We will turn on the 
Planck constant h in section 3. The corresponding equations of motion can be deduced from the action 
with one governing on the electromagnetic field and the other on the scalar field. We write down only 
the second equation of motion for the purpose of deriving the acoustic black hole geometry 

□0 + 2ieA^V^(i) - e^A^A^'cp + - 6 | p = o, (2.2) 

in the d^^A^ = gauge. 

Multiplying (12.21) with 0*, we obtain 

0* [-dtd'ct) + ViVV + 2zeA^V^0 - e^A^A^cf) + 771^00* - 6 | |4= o, (2.3) 

where the index fi = t,xi...X3 and i = X1...X3. We can obtain another equation by multiplying the 
conjugation of (12.21 ) with 

{-d'f(j)* + ViV>* + 2teA^V^'(f)* - e^A^A^cj)*) + 771^00* - 6 | 0. (2.4) 
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The combination of (|2.3I) and (12.41) yields two equations 



- dt {4>*d'4> - 04>*) + Vi VV - 4>V'4>*) + 2ieA^ {^*V^4> + </'V'^(/.*) = 0. (2.5) 



and 



-dt (0*9*0 + 05*0*) + 2|9t0p + V. (0*V> + 0V>*) - 2\di(P\^ + 2m2' 
+2ieA^ (0*V''0 - 0V^0*) - 26101^ - 2eM^A'^|0|2 = 0. (2.6) 



With the assumption = t)e^^^^'*\ the above equations reduce to 



- eA)J + V, [p (V^^ + eA,)] = 0, (2.7) 
e - eAtf - (Vi^ + eA,f + - 6p = 0, (2.8) 



where (12.71) is the continuity equation and (12.81) looks like an equation describing hydrodynamical fluid 
with an extra quantum correction term which is the so called quantum potential. (12.71) and (12.81) are 
actually relativistic and thus are Lorentz invariant. This is a crucial difference from previous literature 
where only the Galilean invariant fluid equations were considered. Hydrodynamics is a long wavelength 
effective theory for physical fluids. Quantum gravity, or even very exotic field theories, when heated up 
to finite temperature, also behave hydrodynamically at long-wavelength and long-time scales. Now we 
consider perturbations around the background (po, Qq)'- P = Po + Pi and 9 = 9q + 9i, and rewrite (12.71) 
and (IZSl) as 



-a 



PoOi + (^0 - e A)pi + V ■ [po V^i + pi ( V^o + eAi)] = 0, (2.9) 



2{9o - eAt)9i - 2{V9o + eAi) ■ V^i - 6pi + = 0, (2.10) 
where we have defined 

1 _3 1 _1 _1 

D2P1 = --Po '(□VA))pi + -Po '□(Po Vi)- 

In this work, we do not consider fluctuations of the electromagnetic field A^. It is a background field that 
can affect the acoustic black hole geometry. The equations (|2.9I) and (|2.10l) corresponds to a relativistic 
fluids equation if bpi term and -D2P1 (the quantum potential) term are absent. Actually, the D2 term 
contains the second derivative of slowly varying p, which can be negligible in the hydrodynamic region 
where k, to are small. We will turn on the quantum potential term only in next subsection. In fact it is 
well-known that in the Madelung representation of the condensation. 
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the Schrodinger equation can be rewritten as a continuity equation plus an Euler equation for an irrota- 
tional and inviscid fluid when the quantum potential can be neglected. The Abelian Higgs model here 
shows that the same property holds for relativistic fluids. 

For the simplicity of calculation, let us introduce variables coq and Vq with the definition 



cuo = —Oq + eAt and vq = V9o + eA. 



We then obtain 



-dt 
-V 



PoV^i 



-V0UJ061 - vo ■ VdiVQ 



0. 



Comparing the above equation with the massless Klein-Gordon equation 

1 



-9 



dAV^g^'d.e,) = 0, 



one may find 



-99 



Defining the local speed of sound 



Po, 



the metric for acoustic black hole yields the form 



9^^u 



bpo 
2cs 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



-v' ': (1 + c2 - t;2)5ii + 

where we made replacement v^coq and notation v'^ = v^v\ We now have a relativistic version 
of acoustic black hole which is a generalization of It is worth noting that the sound velocity Cg is a 
function of the electromagnetic field At, which is treated as a slowly varying quantity. 

In the non-relativistic limit, Cg <^ l,v <^ 1, and keeping only the leading order term of each line 
element, we can recover the metric from non-relativistic theory first obtained by Unruh fill, but with an 
overall factor difference 



9^^u 



(2.16) 
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Equivalently, we can rewrite (|2.15|) as 

bpo 1 



2c. ^/TT^ 
bpo 1 



[-{cl - v'^)dt^ - 2v-dxdt + {v- dxf + (1 + - v'^)dx^] 



2 c, ^ 



-{cl-v^)dr^ + {l + cl-v^) 



V V- 



+ 5ij I dx^dx^ 



(2.17) 



where we have defined a new time coordinate by dr = dt + 



v-dx 



In this coordinate the acoustic 



geometry is actually stationary. 

In the case Vz = , but Vr ^ and v^j, ^ 0, we can simplify the metric in a Kerr-like form 

1 



ds' 



2 c. 



-N {Vr,V^)dl 



-N\Vr,V^ydip^ 



{v^dr — rdipY 



dr"^ + v^l + c2 - v'^dz^ 



where 



0" 



and the coordinate transformations have been used 



dt = dT + 



Vr 



vJ — d 



dr. 



VrV^ dr 



Vr — ct r 



(2.18) 



(2.19) 



(2.20) 



When Vz ^ 0, the metric becomes non-spherically symmetric and non-axisymmetric. In non-relativistic 
limit, the geometry (12.181) reduces to a "draining bathtub" metric first appeared in |[2l. 

At the region far from the acoustic horizon, the normal modes of the 9i field are purely outgoing. 
The operator 9i can be expanded in terms of modes in the region outside the horizon. An observer 
traveling with the fluid as it flows through the acoustic horizon will see that the field 9i is essentially 
the vacuum state in analogy with the behavior of the normal modes of a scalar field in Schwarzschild 
coordinates by a freely falling observer ||T|. The corresponding Hawking temperature at the event 
horizon reads 



h d{cs — f *) 



H ^ r o I horizon- (2.21) 

ZTTKb OXi 

The relativistic version of Hawking temperature from Abelian Higgs model is identical with the non- 
relativistic one. 

Since there is no real gravity, the Newton constant G^r is not involved in the formulation of the 
acoustic "surface gravity" and Hawking temperature. One may compare the acoustic Hawking temper- 
ature with Unruh temperature in Rindler coordinates. The acoustic temperature and the Unruh temper- 
ature yield almost the same form: The acoustic temperature is proportional to the gradients of the fluid 
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kinetic energy (i.e. oc di{cl — Wj^)|horizon ) at the horizon, while the Unruh temperature (T = |^), ob- 
served by an accelerating observer is proportional to the acceleration a^j (= According to Einstein 
equivalence principle, no experiment can distinguish the acceleration due to gravity from the inertial 
acceleration due to a change of velocity. Thus, the Unruh effects in Rindler coordinates and Hawking 
effects in astrophysical black hole spacetime are equivalent. It is surprising that the acoustic geometry 
provides us another way to realize the equivalence principle by the relation (Acceleration) = (Intensity 
of the gravitational field)=(Gradients of the kinetic energy). One point that should be made clear is that 
acoustic analog models of black holes could be used to mimic the kinematical but not the dynamical 
aspects of gravitational systems. For example, it is impossible to reproduce black hole thermodynamics 
using acoustic analogue models of black holes ifTSl . 

2.1 Dispersion relation 

We have shown that by studying the linearized equations of motion of the complex scalar fields of 
Abelian Higgs model, the perturbations obey a Klein-Gordon equation in the curved space-time. Here 
we stress that since we worked in the long-wavelength and slowly varying limit, the quantum potential 
term (the D2 term) is neglible. 

In the high-momentum or eikonal regime, the contribution of the quantum potential would be im- 
portant. In ifTTl . the authors explored the dispersion relation with the contribution of the quantum 
potential. In our case, both the quantum potential and the relativistic effects will contribute to the dis- 
persion relation. IN this regime, the phase fluctuation 9i and the density fluctuation pi can be treated 
as a function with a slowly varying amplitude and rapidly varying phase: 



Ol pi 



9i Pi Re{exp(— ifit 



+ ik-x)} 

01 



Pi 



(2.22) 



(2.23) 



The operator D2 can be approximately written as 



^2Pl 



-2P0 '(°VPo)Pi + 
Pli^f^pof d^pod^pi 



□(Po Vi) 



PiDpo , npi 



2pl 2pl 
□pi 



2p2 2po 



2po 



^Po^(A:2-fiV, 



(2.24) 



7 



where /i = t, X1...X3. We have similar result for 0261. Therefore, we can treat D2 = —^Po {k — ^ ) 
in the eikonal regime. The Klein-Gordon equation (12.91) and (12.101) become 

~ ^PoikkjS'^ - ^'^) - (uJo^ + vikjf = 0. (2.25) 



n = "\ , ^-^k. (2.28) 



2 

With the expression of D2 above, Eg. (12. 251) becomes 

^ - n^) - {ooon + vi^hY + \{k''- = O. (2.26) 

Li the eikonal limit, the last term in the above equation will become dominant and we can neglect the 
other two terms, which gives 

n = k. (2.27) 

On the other hand, in the hydrodynamics regime, we can neglect the D2 term again in the left hand 
of Eg. (12. 251) . To be definite, we choose k^ = k6^^, and define = Vq/uq. Then we easily obtain the 
dispersion relation 

The dispersion relation is quiet different from the result of BH. For large wavelength and in the non- 
relativistic limit, it reads 

n = {cs-v^){l + {cs + v^)/2)k. (2.29) 

The sound velocity in hydrodynamic limit has a correction to the familiar value Cg. The main correction 
is simply due to the fluid motion and the next term is higher order. Since f \ the generalized physical 
velocity in the presence of the electro-magnetic fields, depends on background electric and magnetic 
fields, as one can see from 

= (dJo + eA,)/{-e + eAt), (2.30) 

so is the sound velocity. 

We stress that we assume that the back ground fields change very little within a wave length. If cj. 
and f ^ can be neglected, compared with the speed of light, then the relation return to a linear dispersion 
relation ^ cAk\. 



3 Non-relativistic limit: Acoustic black holes in superconductors 

In the non-relativistic limit, the Abelian Higgs model reduces to the Ginzburg-Landau theory 

2ft|-7A(f,t) = --^S2^(f,t) + a(r)^(f,t) + 6(r) |^(f,t) pV^(f,t), (3.1) 
at 2m* 



where T) = V + m* is the mass of a cooper pair, a{T) and 6(T) are two parameters that depend 
on temperature. We have turned on the Planck constant h. The Ginzburg-Landau coherence length is 
defined by 



Assuming ijjif, t) = p{f,t)e^^^^'^\ we have 

^ + V.(p^)=0, (3.3) 

h— = ^ -a(T)- b(T)p. (3.4) 

at 2m* ^ 2 V ; V ;p 

where we have used the gauge V ■ A = 0, and the velocity is denoted by t? = + ^A. The first term 
in the right hand of (13.41) corresponds to the quantum potential. The above two equations are completely 
equivalent to the hydrodynamical equations for irrotational and inviscid fluid apart from the quantum 
potential. In the long-wavelength approximation, the contribution coming from the linearization of the 
quantum potential can be neglected. 

Linearizing Eqs. (l3.3l) and (13.41) around the background (po, do), we get 

^ + V-(pi^;o + poVei) = 0, (3.5) 

^ + vo-V9, + ^-^p, = 0. (3.6) 
ot m 

Eliminating pi, the above equations lead to 

' -t7o ■ V^i + V ■ poV^i - ttt^vq—^ - — -^iTo ■ V^i = 0, (3.7) 



dt \ h{T) dt h{T) J V b{T) dt h{T) 

where vn = + ^A. One can extract a metric from (13.71 ) because it has the same form as a massless 
Klein-Gordom equation in curved space-time. The non-relativistic version of acoustic metric from (13.71) 
is read as 

ds'^ = J^I - - vl)dt^ - 2vo ■ dm + df- df], (3.8) 
m*Cs I- J 



where Cs = y ^^Po is denoted as the "sound velocity". Actually, one can write the sound velocity in 
terms of the coherent length 

Cs = . (3.9) 
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One should not identify the "sound velocity" defined here with the sound velocity produced by the 
crystal lattices. Here Cg is a function of the density po of superconducting electrons and temperature- 
dependent parameter b{T). 

We would like to do some calculations on the necessary conditions for the formation of an acoustic 
event horizon in Ginzburg-Landau theory. We will show that it would be a tough task to realize acoustic 
black holes by using type I and type II superconductors. The definition of acoustic black hole requires 
^ > 1 in some region. From the London equation 

V'B-^B = 0, (3.10) 

and the Maxwell equation 

VxB = i2oi, (3.11) 

where A(T) = (^^;^^)^^^ and jg denote the penetration depth and the supercurrent, respectively. We 
have a finite magnetic field , say B = (0, By(x), 0), outside a superconducting region. After simple 
computation, we find the supercurrent yields |fT9l 

1 =^ (3 12) 

The critical current density above is only valid for clean type I superconductors. On the other hand, the 
current density has the form 

3sz = ensVos, (3.13) 

where Us = 2\ip'^\ = 2po denotes the superfluid density, and vqs is the velocity of the Cooper pairs. We 
then have 

^0. V2kH{T) H{T) 



(3.14) 



where k = , and Hc2 = 4^ewV(r) ^ V^f^Hc- In the standard textbook, k denotes the G-L param- 
eter, which describes the difference between type I and type II superconductors: < is for type I, 
and K > is for type II. Hc2 denotes the critical magnetic field, and He the thermodynamic critical 
field. From (13.141) . we can see that for type I superconductors, when vq^ = Cg, the superconducting 
phase is broken and return to Normal states. Now we are going to solve (|3.10l) in the high k limit, 
A ^ ^, which is valid for the cooper-oxide superconductors. We will show that even near the normal 
core of the vortex, it is hard for the speed of superconducting electrons to exceed the "sound velocity"- 
Cs- For simplicity, the vortex is assumed to be infinitely long and axially symmetric so that there are no 
z or angular dependence of its field distribution. We then write Eq. (l3.10l) in cylindrical coordinates 

-9,(r9,B) - ^B = 0. (3.15) 
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This solution is given by [[T9] 



B{r) = ^Ko(r/A), (3.16) 

where Kq^t/X) is the zeroth-order modified Bessel function. The current density is derived by the 
Maxwell equation (13.1 II) . 

We are interested in the asymptotic behaviors at small radial distances r ^ A, where Ki[r/X) ^ ^. 
The speed of Cooper pairs is then given by 

^0. = ^ (3.18) 

Use 

We then find 

^ = ^. (3.19) 

Cs r 

The formation of acoustic black holes requires Vqs/cs > 1, i.e. r < -\/2^. But at the core region of 
the vortex r < ^, the order parameter and the superconducting current would decay as r — 0, since it 
corresponds to a normal state. (13.131) tells us that for fixed electron velocity vqs, the current decreases 
as ns{= 2po) decreases. Since the "sound velocity" Cg depends on the density pq, we expect that in the 
region ^ < r < electron velocity can exceed the sound velocity. However, it would be interesting 
if one can verify this point experimentally. 

In order to trap phonons, one need input a sink so that electrons can propagate along 0- and r- 
direction simultaneously. For this purpose, we need turn to spiral flux vortices in the type II supercon- 
ductors. The spiral flux vortices [}20|| in current-carrying type-II superconducting cylinders subjected 
to longitudinal magnetic fields can be used to mimic acoustic black holes. Consider a type-II super- 
conducting cylinder of radius a containing a spiral vortex, in the presence of a longitudinal applied 
magnetic field Ha and transport current /, where Ha and / are positive in the z direction. The Lorentz 
force imposed on the vortex were found to be radially inward. Thus the spiral vortex has a tendency to 
move to the center of the cylinder. One can determine the terminal speed Vq^ of the contracting vortex 
by balancing the Lorentz force and the viscous drag force i^Vq^, where 77 is a phenomeno logical viscous 
drag coefficient per unit length of vortex [[20ll . 

Actually, the spiral vortex in the presence of a sufficiently large current density parallel to its axis, 
is unstable against a helical perturbation ll20l . But this does not prevent the formation of acoustic 
black hole geometry. In the next section, we investigate the spiral vortex geometry from classical 
hydrodynamics in the non-relativistic frame and give the exact metric for a spiral- vortex geometry. 



'The order parameter also reaches its maximum at r = \/2^ lfT9l . 
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3.1 Spiral- Vortex geometry 

In [|2l|, the author constructed a "draining bathtub" black hole metric by a (2 + 1) dimensional Newtonian 
flow with a sink at the origin. An stable vortex solution to the Navier-Stokes equation requires each 
component of the fluid velocity to be non- vanishing. Such an exact solution to Navier-Stokes equation 
was first found by Sullivan [[2T| and later observed in experiments [[22l|. Now, we utilize the Sullivan 
vortex solution of Navier-Stokes equation to model a 4-dimensional acoustic black hole metric, which 
could be realized in experiments. The non-relativistic Navier-Stokes equation is given by 

dtv + v-Vv = -Vp/p + v\7'^v + f/p, (3.20) 

where v is the fluid velocity, p the fluid pressure, p the mass density, v the shear viscosity and / an 
externally specified forcing function. The continuity equation generally accompanies the Navier-Stokes 
equation 

dtp + V ■{pv) = Q. (3.21) 
By defining v = Vip, the Navier-Stokes equation can be rewritten as 



dti^ + -iV^ljy+ ^-z/V'^ = 0, (3.22) 
^ Jo P 

where we have chosen vanishing external force / = 0. Now linearize these equations around a fixed 
background {po^Poii^o) and set p = Po + pi, p = Po + Pi, and ijj = i/jq + "^i- After linearizing the 
equation of motion and the continuity equation, the equations for pi and tpi yield 

^ + V ■ {pivo + poV^i) = 0, (3.23) 

^ + v,.VtPi + ^- uV^tPi = 0. (3.24) 
dt Po 

One can further assume the fluid is barotropic, which implies that p is a function of p only %. The 
barotropic assumption gives 

dp 



Pi = Q^Pi, (3.25) 



After rearranging (13.241) . we have 



Pi = -^Po {dtiJi + vo ■ V^Ai - z/V'V^i) , (3.26) 



^This assumption was first used implicitly in H] (see also 121). 
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Substituting the above linearized equation of motion back into (|3.23l) . we obtain 



-dt 
+V 



dp 



0. 



We can rewrite the above equation in a Klein-Gordon like form 

dp 



dp 



0-21) 



(3.28) 



where the inverse metric is given by 



PoCs 



VnU 



0^07 



(3.29) 



and the local speed of sound is defined by c~'^ = ^. The acoustic metric gf^^, is given in (12.161) . 

What is different from before is that we have third order derivative term of ipi in the right hand 
of (13.281 ). The similar situation was discussed in ^ and Lorentz symmetry breaking was observed. 
Clearly, by writing 'tpi as tpi ~ exp(— icut + ik ■ ic), one see that in the eikonal approximation, the 
right hand of (13.281) will be dominant, while in the hydrodynamic limit, the right hand of (13.281) can be 
ignored. In the following, we will work in long- wavelength approximation and neglect those terms. 

In the cylindrical coordinate (r, 



z), the exact solution to Navier-Stokes equation yields Bill 

2/3 z/ 



—ar + 



r 



_ ar 

e 



To H{x) 
2'Kr Hi oo) 



2az 



l-(3e 



ar 
2u 



where 



flT 



exp 



/3 



dr I ds. 



(3.30) 
(3.31) 
(3.32) 

(3.33) 



To is a constant that depends on the boundary condition, u denotes the shear viscosity, a > and (3 > 1 
are free parameters that can be fixed in experiments. By using the coordinate transformation 



dt = dr 



dip + 



,r2 



-dr 



+ ^,-2 _ o 



dz. 



dr 



.r2 



c? r 



+ 



dz 



ci r 



(3.34) 
(3.35) 
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We obtain a non- spherically symmetric metric by rewriting (|2.16|) as 

jdrdz - 2v^rdipdT + r'^dip'^ } . (3.36) 



It is worth noting that when = 0, we can recover the " draining bathtub" metric obtained in [|2l, that 
is to say 

2 

ds^ = P^\ - (d - vl)dT^ + , . dr"" + dz^ - 2v^rdipdT + r^d^^]. (3.37) 

The form of the metric obtained in (13.361) is non-spherically symmetric and even non axi- symmetric. 
There exists an event horizon for the metric. The event horizon can be fixed by using the Null surface 
equation [[24| 

df df 

a'^jnijr^ = 0, (3.38) 

where f = /{t, r, Lp, z) = denotes the hypersurface. Since the radius of the event horizon would vary 
as a function of r and z, we find that the event horizon is determined by 



where the dot and the prime denote the derivative with respect to r and z, respectively. 



c;^ (3.39) 



4 Conclusion 

In summary, a new acoustic black hole metric has been derived from the Abelian Higgs model. In 
the non-relativistic limit, the metric can reduce to an acoustic black hole geometry given by Unruh. 
Since the Abelian Higgs model describes high energy physics, our results indicate that acoustic black 
holes might be created in high energy physical process, such as quark matters and neutron stars. The 
dispersion relation was given and in the non-relativistic limit, we can recover the linear dispersion 
relation. We have interpreted the requirements of creating an acoustic black hole by using the language 
of superconductor physics. Our results have demonstrated that although it would be very difficult 
to mimic the acoustic black hole in type I superconductors, people would be able to to discover the 
acoustic black holes near the magnetic vortex in type II superconductors . We developed an acoustic 
black hole with spiral vortex geometry from the Navier-Stokes equation in classical Newtonian fluids. 

One point to be figured out is that we have used the vortices to mimic the black hole. Actually, 
when a superfluid flow exceeds the Landau critical velocity, vortices can be created and excitations 
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with negative energy appear, which means that black hole configurations are energetically unstable 
since the speed of sound is also the the Landau critical velocity. However, this does not mean that 
the supercritical flow of the superfluid vacuum is not possible. In superfluid ^He the Landau velocity 
for vortex nucleation can be exceeded by several orders of magnitude, without vortices creation. In 
[|23l . the event horizon was proposed to be realized in superfluid ^He by arranging the superfluid flow 
velocity to be constant but the "sound velocity" to change in space and in some region becomes less 
than the flow velocity of the liquid. In our above discussions, we have shown that in the ideal type 
II superconductors, the "sound velocity" (13.91) indeed depends on the superfluid density po, but the 
superconducting (superfluid) flow velocity can be chosen to be constant by varying the current and 
the superfluid density to be po simultaneously in (13.131) . Therefore, if the "sound velocity" changes in 
space, superconducting (superfluid) flow velocity would exceed the "sound velocity" in some regions. 
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